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Abstract 

The recoil correction to the ground state energy of hydrogenlike 
atoms is calculated to all orders in aZ in the range Z = 1-110. The 
nuclear size corrections to the recoil effect are partially taken into 
account. In the case of hydrogen, the relativistic recoil correction 
beyond the Salpeter contribution and the nonrelativistic nuclear size 
correction to the recoil effect, amounts to -7.2(2) kHz. The total 
recoil correction to the ground state energy in hydrogenlike uranium 
(238tj91+) constitutes o.46 eV. 
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1 Introduction 

The complete aZ-dependence formulas for the nuclear recoil corrections to 
the energy levels of hydrogenlike atoms in the case of a point nucleus were 
first obtained by a quasipotential method [1] and subsequently rederived by 
different approaches [2-4]. According to [4], the nuclear size corrections to the 
recoil effect can be partially included in these formulas by a replacement of 
the pure Coulomb potential with the potential of an extended nucleus. The 
total recoil correction for a state a of a hydrogenlike atom is conveniently 
written as the sum of a low-order term AEl and a higher-order term AEjj 
[1], where (h — c — 1) 

A£l = ^(a|[p 2 -(D(0).p + p.D(0))]|a), (1) 

AE H = -J— f°° du{a\(B(u)-^h 

xG(u + e a )(B(u) + ^h\a). (2) 

Here, \a) is the unperturbed state of the Dirac electron in the nuclear po- 
tential V(r), p is the momentum operator, G(uo) = (ou — H(l — i0))~ l is the 
relativistic Coulomb Green function, H — (a ■ p) + f3m + V , a\ (I — 1, 2, 3), 
(3 are the Dirac matrices, 

D m (uj) = -4naZaiDi m (u) , 

and Dik{u), r) is the transverse part of the photon propagator in the Coulomb 
gauge: 

, . 1 rexp (i\uj\r) (exp (i\uj\r) — 1) •. 

47r L r u z r > 

In equation (2), the scalar product is implicit. For point-like nuclei, V(r) = 
Vcij') = —ctZ/r. If extended nuclei are considered, V(r) is the potential of 
the extended nucleus in eq. (2) and in calculating e a , \a), and G{uj). There- 
fore, the nuclear size corrections are completely included in the Coulomb 
part of the recoil effect. In the one-transverse-photon part and the two- 
transverse-photon part (see Ref. [4]), they are only partially included. At 
least for high Z we expect that this procedure accounts for the dominant 
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part of the nuclear size effect since using the extended nucleus wave function 
and the extended nucleus Green function strongly reduces the singularities 
of the integrands in (1) and (2) in the nuclear region. 

The term AEl contains all the recoil corrections within the (aZ) A m 2 /M 
approximation. Its calculation for a point nucleus, based on the virial rela- 
tions for the Dirac equation [5-7], yields [1] 

AE L = m2 ~ £ "° , (3) 
L 2M ' K J 

where e a o is the Dirac electron energy for the point nucleus case. 

AEh contains the contribution of order (aZ) 5 m 2 /M and all contributions 
of higher order in aZ. To lowest order in «Z, this term represents the 
Salpeter correction [8]. The calculation of this term to all orders in aZ 
was performed in [9,10] for the case of a point nucleus. According to these 
calculations, the recoil correction to the Lamb shift of the Is state in hydrogen 
constitutes -7.1(9) kHz, in addition to the Salpeter term. This value is close 
to the (aZ) & m 2 / M correction (-7.4 kHz) found in [3] and is clearly distinct 
from a recent result for the (aZ) 6 m 2 /M correction (-16.4 kHz) obtained in 
[11]. (The (aZf log (aZ)m 2 /M corrections cancel each other [12,13].) The 
total recoil correction to the ground state energy in 238 U 91+ was calculated 
in [9] to be 0.51 eV. 

In this work we calculate the recoil correction to the ground state energy 
of hydrogenlike atoms in the range Z— 1-110 using the formulas (1) and (2) 
employing the potential of an extended nucleus. 



2 Low-order term 

Using the virial relations for the Dirac equation in a central field [7], the 
formula (1) can be transformed to (see Appendix) 

+2aZK(a\a z 5V/r\a) - 2e a (a\5V\a) 

+2(m + 2e a K,)(a\a z 5V\a) — 2aZm(a\a x SV\a) 

-Ame a {a\a x r5V\a)} , (4) 
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where e a and e a o are the Dirac electron energies for an extended nucleus and 
the point nucleus, respectively, k = (— l) J+ ' +1 / 2 (j + 1/2) is the relativistic 
angular quantum number, 5V = V(r) — Vc(r) is the deviation of the nuclear 
potential from the pure Coulomb potential, and a x and a z are the Pauli 
matrices. Here, the notations for the radial matrix elements from [7] are 
used: 

(a\u\b) = / [G a (r)G b (r) + F a (r)F b (r)]u(r)dr, (5) 
Jo 

roo 

(a\a z u\b) = / [G a (r)G b (r)-F a (r)F b (r)]u(r)dr, (6) 
Jo 

roo 

(a\a x u\b) = / [G a (r)F b (r) + F a (r)G b (r)}u(r)dr. (7) 
Jo 

G/r — g and F/r = f are the radial components of the Dirac wave function 
for the extended nucleus, which are defined by 



,/, ( r \ — ( 5 , nfc( r )^Km( n y 

^W(r) - ^ tfnfi{r)Q _ Km{n) 



The first term on the right side of equation (4) corresponds to the low-order 
recoil correction for the point nucleus (see Eq. (3)). The second term gives 
the nuclear size correction . We calculate this term for the uniformly charged 
nucleus. In Table I, we display the results of this calculation for the the Is 
state. The values are expressed in terms of the function AFl(o:Z) which is 
defined by 

AE L = ^^(l + AF L (aZ)). (8) 

In order to compare the nuclear size correction to the low-order term with 
the corresponding correction to the higher-order term (see the next section), 
in the last column of the Table I we display the value AP L (aZ) which is 
defined by 

Using Eq. (4), one easily finds for an arbitrary ns state and for very low 
Z (aZ < 1) 



1 



n 



12, sQ 72 



AF L (aZ) = -(aZy(Rmy - —{aZfRm , (10) 



5 v 7 v ' 35 
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where R = ^5 /3(r 2 ) 1 ^ 2 is the radius of the uniformly charged nucleus. The 
first term in (10) is a pure nonrelativistic one. It describes the reduced mass 
correction to the nonrelativistic nuclear size effect. So, if the nuclear size 
correction to the energy level is calculated using the reduced mass, this term 
must be omitted in equation (10). The second term, which is dominant, arises 
from the Coulomb part ((a|p 2 |a)/(2M)). For the standard parametrization 
of the proton form factor 



ftp) 



which corresponds to 



p{r) 



and 



A 4 



(A 2 +p 



2\2 ' 



8tt 



(11) 



exp (— Ar) 



(12) 



aZ 

V(r) = 

r 



-exp, 



the contribution of this term to AP L is 



-Ar)(2 + Ar) 



. „, 35 m 
L 8 A 



(13) 



(14) 



We will see in the next section that this term cancels with the corresponding 
correction to the Coulomb part of the higher-order term. This implies that 
the sum of the low-order and higher-order contributions is more regular at 
r — > than each of them separately. 



3 Higher-order term 

To calculate the higher-order term (2) we transform it in the same way as it 
was done in [9]. The final expressions are given by the equations (41)- (54) of 
Ref. [9] where the pure Coulomb potential (Vc(r) = —aZ/r) in the equations 
(42) and (48) has to replaced by the potential of the extended nucleus V(r). 
We calculate these expressions for the uniformly charged nucleus by using the 
finite basis set method with the basis functions constructed from B-splines 



5 



[14] . The algorithm of the numerical procedure is the same as it is described 
in [9]. The results of the calculation for the Is state are presented in the 
second column of the Table II. They are expressed in terms of the function 
P(aZ) defined by 

For comparison, in the third column of this table we list the point-nucleus 
results (Po(aZ)) that are obtained by the corresponding calculation for R — > 
0. These point-nucleus results are in good agreement with our previous 
results from [9]. In the fourth column of the table, the difference AP = P—Pq 
is listed. Finally, in the last column the Salpeter contribution [8,15] 

if S) = -| In (aZ) - 8 - 2.984129 + H l n 2 + f (16) 

is displayed. 

For low Z the nuclear size correction to the higher-order term is mainly 
due to the Coulomb contribution 

AE « ] = 2^/I^FT7oP M<a||p ' V1G(w + £j[p ' K||a> - (17) 

It is comparable with the deviation of the complete aZ-dependence value 
from the Salpeter contribution (in the case of hydrogen AP = —0.0092(2) 
while Pq — Ps = —0.0162(3)). To check this result let us calculate the finite 
nuclear size correction to the Coulomb part of the (aZ) 5 m 2 /M contribution. 
Taken to the lowest order in aZ, formula (17) yields 

,(C) _ (27r) 3 ! ^ y/p 2 + m 2 -m p 2 V 2 (p) 



AE « = "WW) / / 2 ' 2 + V2 / 2. 2 > 

2M J (VP + m + m ) VF + m 

where a (O) is the non-relativistic wave function at r = and V(p) is 
the nuclear potential in the momentum representation. Using the standard 
parametrization of the proton form factor 

-Tr, \ OiZ A 4 

V(p) = - n 9 9 , 9 . 9 19 
2n 2 p 2 (A 2 +p 2 ) 2 
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and separating the point nucleus result from (18), we can write for an ns 
state 



[aZfm 2 
nn 3 ~M 



;-4/3 + AP (C) ) , 



where 
AP (c) 



P 



m 



-4 / dp 

Jo Wp 2 + m2 + m ) 3 Vp 2 + m ' 2 ^ (A 2 + P 2 ) A 



A 8 



- 1 



(20) 



(21) 



Evaluation of this integral to the lowest order in m/A yields 

{c) = 35 m 
8 A 



(22) 



As we noted above, the correction (22) cancels with the corresponding cor- 
rection to the low-order term (see Eq. (14)). For (r 2 ) 1 / 2 = 0.862(12) fm [16], 
which corresponds to A = \f\2/{r 2 ) 1/2 = 0.845 m p = 793 MeV, the formula 
(22) yields AP^ = 0.00886 while the exact calculation of the integral (21) 
amounts to AP^ = 0.00874. These results are in good agreement with the 
corresponding result (AP = 0.0092(2)) from the Table II. 



4 Discussion 

In this work we have calculated the recoil correction to the ground state 
energy of hydrogenlike atoms for extended nuclei in the range Z — 1 — 110. 
This correction is conveniently written in the form 

AE-— + —-—P FS (aZ). (23) 

The function P FS (aZ) = P(aZ) + AP L (aZ) is shown in Fig. 1. For compar- 
ison, the point nucleus function P (aZ) and the Salpeter function P$(aZ) 
are also presented in this figure. The Table III displays the values of the 
recoil corrections (in eV) in the range Z— 10-1 10. 

In the case of hydrogen we find that the recoil correction amounts to 
AE = —7.2(2) kHz beyond the Salpeter contribution and the nonrelativistic 
nuclear size correction to the recoil effect (the first term in Eq. (10)). It 
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almost coincides with the point nucleus result. This is caused by the fact 
that the nuclear size correction to the higher-order term (Eq. (22)) and the 
relativistic nuclear size correction to the low-order term (Eq. (14)) cancel 
each other. 

For high Z, where the aZ expansion as well as the reduced mass approx- 
imation are not valid any more, we should not separate any contributions 
from the total recoil effect. In the case of hydrogenlike uranium ( 238 U 91+ ), 
the total recoil correction constitutes AE = AE L + AE H = 0.46 eV and is 
by 10 % smaller than the corresponding point nucleus value (A_E p . n . = 0.51 
eV) found in [9]. This improvement affects the current numbers of the Lamb 
shift prediction [17]. 

Finally, we note a very significant amount of the nuclear size effect for 
Z—110. According to the Table III, the finite nuclear size modifies the point 
nucleus result by more than 40%. 
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Appendix 



Using the identity p 2 = (a ■ p) 2 , the Coulomb part of the low-order term can 
be written as 

A4 C) = (a\j^\a) = ^ M (a\(e a -(3m-V) 2 \a) 

= ^l + m 2 + (a\{V 2 -2e a V)\a) 

+ 2m{a\P(V-e a )\a)]. (24) 

As described in detail in [18], the Breit part of the low-order term can be 
transformed to 

A4 B) = -^(a|[D(0)-p + p.D(0)]|a> 
1 . .aZ / (a • r)\ 

n (a\-(2e a -2(3m-2V + — a r (3)\a) 



2M X 'rV r 
1 

2M 



— [2aZ(a\V/r\a) - 2aZe a (a\l/r\a} 

+ 2aZ(a\m(3/r\a) + 2naZ \ g a f a dr, (25) 

Jo 

where a r — (a • r)/r and k = (— l) J+ ' +1 / 2 (j + 1/2) is the relativistic angular 
quantum number of the state a. In the following we will use the notations of 
Ref. [7], 

A s = / (G 2 + F 2 )r s dr , (26) 
Jo 

/*oo 

5 s = / (G 2 - F 2 )r s dr , (27) 

JO 

/•oo 

C s = 2 GFr s dr, (28) 

where G/r = (7 and F/r = f are the radial components of the Dirac wave 
function for the extended nucleus, and the radial scalar product defined by 
the equations (5)-(7). Using the equation (2.9) of Ref. [7], we find 

AE L = -±-iel- m i + (a\5V(5V-2e a )\a) 

-aZ(aZA- 2 - kC' 2 - 2mB~ 1 )} , (29) 
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where SV = V - V c = V + aZ/r. /From the equations (2.8)-(2.10) of Ref. 
[7], one obtains 



aZA- 2 - kC- 2 
laZmB- 1 



— 2n{a\a Z 8V / r\a) + 2m(a\a x SV\a) , 
2(m 2 - e 2 a ) + 2(m + 2e a K){a\a z 5V\a) 
— Ame a {a\a x r5V\a) . 



Substituting (30) and (31) into (29), we find 

1 



AE r 



m 2 -e\ 



. y a\(5V) 2 \a) + 2aZK(a\a z 5V/r\a) 
— 2e a (a\5V\a) + 2(m + 2e a K)(a\a z 5V\a) 
—2aZm(a\a x 5V\a) — Ame a (a\a x r5V\a) 



(30) 
(31) 



(32) 



Separating the point nucleus result from the right side of (32), we get the 
equation (4). 
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Table 1: Nuclear size correction to the low-order term for the Is state ex- 
pressed in terms of the functions AFi(aZ) and AP^aZ), defined by equa- 
tions (8) and (9), respectively. The values of the nuclear radii employed in 
the calculation are taken from [16,19-23]. 



z 


(r 2 ) 1 / 2 , fm 


AF L (aZ) 


AP L (aZ) 


1 


0.862 


-0.337xl0- 8 


-0.0136 


2 


1.673 


-0.519xl0~ 7 


-0.0262 


5 


2.397 


-0.102xl0~ 5 


-0.0329 


10 


3.024 


-0.976xl0- 5 


-0.0394 


20 


3.476 


-0.933xl0~ 4 


-0.0472 


30 


3.928 


-0.406xl0~ 3 


-0.0607 


40 


4.270 


-0.126xl0~ 2 


-0.0797 


50 


4.655 


-0.340xl0~ 2 


-0.1099 


60 


4.914 


-0.823xl0- 2 


-0.1539 


70 


5.317 


-0.0195 


-0.2295 


80 


5.467 


-0.0436 


-0.3442 


90 


5.802 


-0.0993 


-0.5506 


92 


5.860 


-0.117 


-0.6073 


100 


5.886 


-0.224 


-0.9038 


110 


5.961 


-0.517 


-1.572 
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Table 2: Higher-order term for the Is state expressed in terms of the function 
P(aZ) defined by Eq. (15). The nuclear radii employed in the calculation 
are the same as in Table I. P (aZ) is the related value for the point nucleus 
and AP = P — P . P s (aZ) is the Salpeter contribution obtained by Eq. 
(16). 



z 


P(aZ) 


Po(aZ) 


AP(aZ) 


P s (aZ) 


1 


5.4391(3) 


5.4299(3) 


0.0092(2) 


5.4461 


2 


4.9703(3) 


4.9528(3) 


0.0175(2) 


4.9840 


5 


4.3281(3) 


4.3034(3) 


0.0247(2) 


4.3731 


10 


3.828 


3.795 


0.031 


3.9110 


20 


3.330 


3.294 


0.036 


3.4489 


30 


3.086 


3.044 


0.043 


3.1786 


40 


2.977 


2.927 


0.050 


2.9868 


50 


2.973 


2.914 


0.060 


2.8380 


60 


3.072 


3.006 


0.066 


2.7165 


70 


3.295 


3.234 


0.061 


2.6137 


80 


3.686 


3.672 


0.013 j 


2.5247 


90 


4.330 


4.521 


-0.191 


2.4462 


92 


4.501 


4.779 


-0.277 


2.4315 


100 


5.40 


6.41 


-1.01 


2.3759 


110 


7.24 


12.43 


-5.19 


2.3124 
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Table 3: Recoil corrections in eV. For comparison, the nonrelativistic recoil 
correction is given separately. The last column displays the deviation from 
the point nucleus results for the total recoil effect. The mass values are given 
in nuclear mass units. They were taken from [24], except for Z — 110 where 
we adopted the value of [21]. 



z 


M/A 


nonrel. recoil 


total recoil 


finite size effect 


10 


20.2 


0.037 


0.037 




20 


40.1 


0.075 


0.075 




30 


65.4 


0.104 


0.105 




40 


91.2 


0.134 


0.137 




50 


118.7 


0.163 


0.171 




60 


144.2 


0.196 


0.215 


-0.001 


70 


173.0 


0.227 


0.269 


-0.003 


79 


197.0 


0.26 


0.33 


-0.01 


80 


200.6 


0.26 


0.34 


-0.01 


82 


207.2 


0.27 


0.36 


-0.01 


90 


232.0 


0.30 


0.44 


-0.03 


92 


238.0 


0.30 


0.46 


-0.05 


100 


257.1 


0.34 


0.61 


-0.14 


110 


268.0 


0.42 


0.97 


-0.75 
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nuclear charge number Z 

Figure 1: The function P FS (aZ), compared to P (aZ) and P s (aZ). 



